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Abstract 

A Green's function formalism to analyze the scattering properties in confined geometries is de- 
veloped. This includes scattering from a central field inside the guide created e.g. by impurities. 
For atòmic col·lisions our approach applies to the case of parabòlic confinement and, with certain re- 
strictions, also to an anharmonic one. The coupling between the angular momentum phase shifts 5i 
of a spherically symmetric scattering potential V(r) due to the cylindrical confinement is analysed. 
Under these general conditions, a broad range of scattering energies covering many transversal ex- 
citations is considered and changes to the bound states of V(r) are derived. For collisions between 
identical àtoms, the boson-fermion and fermion-boson mappings are demonstrated. 

PACS numbers: 03.75.Be, 05.30.Fk, 05.30.Jp, 34.10,+x 
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Implementing atom-optical devices often requires a strong confinement for all except one 
degree of freedom |1, and Refs therein]. Examples of physical situations where a strong 
confinement is needed are guided matter-wave interferometers lli, one dimensional optical 

_i- „ a — 

Girardeau gas |4J or phase fluctuations of quasi-condensates p| . A proper description of the 
dynamics of such reduced quasi-lD systems should account for the nature of the discrete 
transverse states. Therefore one needs to deduce the effective 1D interaction between the 
remaining longitudinal degrees of freedom from the real 3D free-space interactionpotential 
V{r). Col·lisions under confinement different from the 1D case are treated in jÜQ]. Apart 
from ultracold atom-atom col·lisions, scattering in confined geometries also occurs in various 
physical situations such as scattering of guided atòmic matter waves or of guided electro- 
magnètic and accoustic waves [8( from obstacles inside a guide, e.g., (heavy) impurity àtoms 
or material defects, respectively. The latter is of importance for the propagation of radiation 
or sound within transmission lines or resonators. 

As for atom-atom col·lisions, resonant quasi- 1D scattering in the transverse ground state 
of the guide (single mode regime) was first considered for bosons in an harmònic guide 
employing for the interaction potential V(r) a delta-like zero-range approximation 9]. Nu- 
merical simulations 10] confirmed for certain finite range potentials V(r) the existence of the 
so-called confinement induced resonance (CIR) originally predicted in 9|. A further investi- 
gation of the CIR is provided in [Uj dealing for the first time with a general finite-range V(r) 
for both bosons and fermions under harmònic confinement. Effects of the non-parabolicity 

n 

of the confinement are considered in 1 2J , with focus on the center of mass dynamics and 
employing a zero-range approximation for the interaction V(r) . 

The present work extends the above approaches and gives an alternative and complemen- 
tary description of scattering under confinement, treating both the cases of col·lisions and of 
scattering by fixed obstacles. We develop a general formalism based on the Green's functions 
that allows us to express the scattering properties in confined geometries in terms of the 
phase-shifts à~i of free-space scattering. The coupling between these phase-shifts is explicitly 
taken into account. A general initial scattering state can be treated properly, describing 
in particular the "multi-channel" regime, in the sense that the total energy allows several 
transversal excited states to be effectively occupied. 

In the case of col·lisions where V(r) is the atom-atom interaction potential, the center of 



2 



mass motion is known to separate from the relative one only for a parabòlic confinement. 
Our approach then provides a deeper understanding of this collision process. On the other 
hand, for an arbitrary confinement, scattering processes triat can naturally be described 
by the formalism include, e.g., the quantum scattering of individual cold àtoms, or other 
equivalent systems, by a central field V(r) fixed in the center of the guide at r — 0. As for 
atom-atom scattering, the relative coordinates r are not exactly separable from the center 
of mass coordinates R if the confinement is no longer parabòlic. Nevertheless, in such 
a situation of coupled center of mass and relative motion, the formalism provides in the 
ultracold regime a distinct starting-point to account for this coupling. 

Under the above restrictions concerning atom-atom col·lisions, our investigation confirms 



that the CIR 01 is a general consequence of the dominant terms of the scattering amplitudes. 
The main requirements are a large positive s-wave scattering length a, a ~ /j_ is the 
length scale of the confining potential U(p) , such that U ~ for p <C l±, and equals 
the cylinder radius for a square-well type confinement], a short-ranged scattering potential 
V(r), Ry <C l± [Rv is the range of V(r) , such that V ~ for r ^> Ry], small longitudinal 
momenta and small phase-shifts ó~i, as described below. The resonance is accompanied by a 
l = bound-state of V(r) strongly distorte d by the confinement U (p) and pushed towards 
the continuum. This modified bound state [l0| is shown to be a herald of the CIR. In the 
context of scattering of individual guided àtoms by a central field, these conclusions hold 
irrespective of restrictions due to anharmonicities and imply the unambiguous strong effects 
of confinement on the scattering process. 

Phase-Shifts. The Schròdinger equation for the scattering wave function ^(r), with 
r = (p, z), reads 



where u(p) = 2pU(p)/H 2 , v(r) = 2pV(r)/h 2 , and E = h 2 k 2 /2p > is the total energy. 
In the case of atòmic col·lisions, r is the relative coordinate and the relation between U(p) 
and the confining potential U c {pi) of the z-th particle in the laboratory reference frame is 
given by U(p) = 2U c (p/2) . Note that this relation is no longer exact for non-parabolic U c 
(but provides a first uncoupled description of the relative motion, by quenching the center 
of mass at the origin R = 0). The cylindrical boundary condition is met by expanding the 
solution in the transverse eigenstates ip n (p), with energies e n = h 2 q 2 /2p and normalized to 




(1) 
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J dxdy (p n (p)*ip m (p) = 5 nm . As a result, one obtains the integral equation 

*(r) = * 4 (r) - / d 3 r'G c (r,r')v(r')y(r'). (2) 



For a given low enough such that k ~ let tie be the integer obeying fc 2 = g 2 B + fe 2 B < 
q\ +nE - The following study includes the situations of ground state scattering (ue — 0) as 
well as scattering in the transversally excited modes (n^ > 1). In both cases, transverse 
states with n > % can only be virtually occupied, since k 2 < q 2 . The general initial state 
is #i(r) = Y^n=o^ lknZ Vni.p) for some constants b n , with A; 2 = g 2 + A; 2 . In Eq.(J2J, 

oo 

G c (r, r') = Vn(p)MpTGn(z - z') (3) 

n=0 

is an axially symmetric Green's function and G n (z) = — e íkn ^/2ik n (for n < Ue) and 
G n (z) = e~ Pn ^/2p n (with k 2 = q 2 — p 2 , for n > 1 + n E ) are 1D Green's functions. The 
excited states with quantum numbers larger than n E decrease exponentially with increasing 
distance from the scattering region. In the asymptotic limit \z\ — > oo, one has for n < n E 

n E 

*(r) « [b n e lKz + /± e lk ^} <p n (p), z -+ ±oo , (4a) 

n=0 

/n± = À~J dr ' [ e± ^*WpO] VW), (4b) 

where is the n-th channel effective 1D scattering amplitude for forward z > and back- 
ward z < scattering. 

Consider next G c (r, r') for r' < r <C In this region t/(p) ~ and one should be able 
to approximate G c by the free 3D Green's functions Gi 2(7", r') = e ±ïfc ' r_r '/47r|r — r'\. Thus, 
we write 

1 r / „ik\r— r'\ p —ik\v—r'\ 

GJr,r') = — / d<j/ 7+ -A r + 

27T 7 V 47r|r — r'| 47r|r — r' 

+ A c (r,r') (5a) 
k^^ji(kr') / ~f + h\ 1 \kr) — 7_/i| 2 ' ) (fcr) 

2/ -I- 1 

-P,(cos0)Pj(cos0') 



/ 
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+ A c (r, r') , r' < r < Z ± . (5b) 

In Eq.()5b|). we have used the well known expansion of Gi^ in spherical coordinates 
Prob.7.5]. The value of j± and A c can be explicitly obtained if U(p) is approximated by a 



square-well type confinement for r <C l±. Indeed, the eigenstates are then close to Bessel 
functions, (f n {p) ~ N n J (q n p) / V 1 / 2 ^, normalized on a disc of radius l±, N n = l/\Ji(r n+1 )\, 
r n+ i being the (n + l)-th root of Jo- Separating from G c the terms n < ng, the series 
for n > he can be approximated by an integral over q emerging from the continuum limit 
q n — ► q and vàlid when r',r <C l±. Note that q starts at qi+n E > k. One then compares 
real and imagjnary parts of G c in Eq(jïïJ) with a suitable expansion of G12 in cylindrical 
coordinates [13, Prob.7.9] in Eq([üa|. This comparison leads to 

n E 

7± = l/2± 7 /2, 7 = ^2iy n 2 /^ n /i, (6a) 

n=0 

A c (r,r') = -±-J PC dpe-^-*y (qp)J (qp f ) , (6b) 

with g = a/A; 2 + p 2 and gi+ nE = \/k 2 + p 2 . The homogeneous (Helmholtz) term A c corrects 
the Green's function 7+Gi +7_G 2 , with 7+ + 7_ = 1, in order to account for the discreteness 
due to the confinement. Within the flatness condition, the above approach is vàlid for 
arbitrary U(p). It yields an intrinsic connection between the confined and the free space 
scattering approaches (see [9| for parabòlic confinement). 

In order to obtain the scattering phases Si that are associated with the spherical symmetry, 
we expand the incident state in spherical coordinates employing e tknZ ip n (p) = J^i^fà + 
í)&niji(kr)Pi(cos9), with a n i = N n Pi(k n / k) /n 1 ^ 2 ^ jï^j]. Analogously in A c , the equivalent 
expansion is given by e~ pz Jo{qp) = Ylt^fàl + l)Pi(ip/k)ji(kr)Pi(cos6) stemming from an 
analytic continuation into the complex #-plane (8 — > n/2 — i9). Inserting these expressions 
and Eq.([5b|) into Eq.(J2J) and using Eq. íjüaf) yields, for R v <C r <C l±, 

*(r) w J2 il ( 2l + 1 )l a i + ^( z )-^ kT ^^ kr ) P i( cosd ) 
1 

+ ^(2/ + !) [kTi]ni(kr)Pi(cos9) , (7) 
1 

with ai = Y2ío b nam. Here W { = i~ l J d z r'[ji{kr')Pi{cosB')}v{r')^{r') and 4ttji{z) = 
j£ c dp j {z) d 3 r'Pi{±ip/k)e ±pz ' J (qp')v(r')^(r'). The integration over r' for 7/(2) is performed 
in a finite volume íl covering the range of v(r'). If z is outside íl, the positive sign refers to 
a positive z and vice-versa. Inside íl, both signs are needed according to whether z ^ z'. 
Except for this z-dependence of ji(z) in Eq.(JZJ), we have now succeeded in representing the 
total scattering wave function in spherical coordinates. 



5 



Noteworthy at this point is the fact that 7z (z) accounts for couplings between different 
angular momenta. Indeed, by using e ±pz> J (qp') = J2i> i 1 ' (21' + l)Pi'(Tip/k)ji'(kr')Pi> (cos 6') 
and the property P//(=fm) = (—) l 'Pi>(±u) in the definition of 711(2), one gets a constant 
7/(2) if, for each /, only Z'-waves are kept such that l + V = even. The latter condition is 
also necessary to obtain non-zero matrix elements (l \ U(p) \ V) due to the parity symmetry 
r — > —r. Therefore, a constant ji(z) ~ 7; arises 

7z = £}(2Z' + l)P lv T v , Z = 0,1,2,... , (8) 

í'M 

where P«/ = k J^ c ^ k du Pi(iu)Pi'(iu) and Z'[Z] denotes the sum over even (odd) /' for even 
(odd) Z. Eq.(jHJ) is equivalent to the condition that the "perturbation" U(p) to the free space 
scattering does not couple even and odd angular momenta. 

It is now possible to introduce the phase-shifts 5i. The solution Eq.(J7J) can be written as 
(R v < r < l±) 

\Er(r) « [cosí; jí(/cr) — siní/ n;(fcr)] P/(cos0), (9a) 

c , = (2l + l)(a l + n/i)i l ^ = q f + 7; 
' cos — Í7 sin 5; ' fyA; — cot 5/ ' 

where the last two relations define formally c[ and 5i. That this 5i is the actual phase-shift 
can be seen as follows. On one hand, Eq. (jüajl is the (intermediate) asymptotics Ry <C r <C l± 
of the solution \&(r) = £\ cJi?/(r)P^ (cos 0) in the region of V(r). On the other hand, the 
free-space scattering solution in this region, i.e., not taking into account the boundary, is 
just a different superposítion ty 3 £)(r) = qPqPj with the same radial part Ri. In other 
words, the effect of the confinement U (p) is to change the superposition coefficients from q 
to c[ while keeping the scattering phases of the free-scattering problem. Then the second 
relation in Eq. (|9bJ) together with Eq.fjHJ) gives a matrix equation for T\ in terms of ói, i.e., 
for 1 = 0,1,2,... 

{ijk - k cot 61) T t = a t + J2( 21 ' + l )Pi' T v ■ ( 10a ) 

/'[/] 

Finally, the effective amplitude f^ is given by expanding e ±lknZ ' <p n (p') m the integrand of 
Eo. pSjl . thus 

fn = fng i /nu = ( ^ ] i ^ ] 

V even i odd 



(« + 1)4^ (10b) 
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The relationship between the amplitudes in Eq. (jlObj) and the matrix elements T} of Eq. ()10a|) 
constitutes the main result of our formalism. 

Current Conservation. Inserting Eqs. ()l()ajl()b|) into Eq.(0J), the probability conservation 
should follow. From the total current along the z-axis, the conservation condition is 

n E 

= £(|/n fl | 2 + MKfng} + \fnu\ 2 + Re{b*J nu })k n . (10c) 
n=0 

In the remainder of this paper, we analyse the scattering process given by the leading terms 
of Eqs. fjlOj) . We consider first the case of the single mode regime in more detail, followed by 
the case of transverse excitations and angular momenta couplings. 

Single Mode Resonances. When only the ground state {n E = 0, b n = 6 0n , k 2 = q$ + k^) 
represents an open channel, the symmetric and antisymmetric sectors of Eq. ()4aj) . ^/(r) = 
[Í> g {z) + ipu(z)} fo(p), are given respectively by (for z ^ 0) 

ifí g (z) = (1 + f 0g ) cos (k z) + i f 0g sia(k \z\), (lla) 
if) u (z) = i(l + f ou) sm (k z) ± f 0u cos (k z). (11b) 

In the context of collisions between identical partides, it is clearly seen that, at resonance 
fo g = —1, the bosonic sector ip g is mapped into a non-interacting f 0u = pair of (spin- 
polarized) fermions, the well known fermionization of impenetrable bosons. Now, the inverse 
is also seen to occur for ip u at the fermionic resonance, fo u = —1, first obtained in |llj]. A 
further insight is gained by setting 

fog,v = ~ t 1 + icotóg^r 1 ■ (12) 

The conservation condition Eq. íjlOcj) is then fulfilled for real 1D phase-shifts 5g )U and one 
can rewrite ip g = e l5g cos (k \z\ + S g ) and ip u = ie tSu sin (k z ± 5 U )- Thus at resonance \S 9iU \ = 
7r/2 and the above discussed boson-fermion and fermion-boson mappings exist also under 
longitudinal confinement, e.g., by imposing Tp gtU (z = In) = 0, as numerically verified in 
Ref. ÏJ. 

CIR and bound-states. The resonance fo g = — 1 can be calculated from a general potential 
V(r) by solving Eq. ()l()a|l for even /. Since kRy y Rv/l± 1, the phase-shifts tan ói = 
tan 5i(k) ~ k 2l+1 ~ l//^ +1 are generally small 14 1 for large l±. From Eq. ()lüa|) . it follows 
that l = is the leading contribution and fo g has the form compatible with Eq. (jl2|) 

fo « fog » 1— » ( 13a ) 



l + l 



ko 
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where d± = l±/N , P 00 = p c , and a is the 3D s-wave scattering length, kcotó ~ —1/a. 
This corresponds to solving for z under an effective 1D pseudopotential Vid(z) = gioS(z), 
with the coupling strength 



Ü 2 2a f Ca 
p d\ \ d 



9w = — -w ( 1 ~ ~T~ ) ' C ' = d ±Pc ■ (13b) 



As in previous works in the single mode regime (for atom-atom col·lisions in parabòlic con- 
finement) , the resonance \gm\ ~ > oo at d± ~ Ca requires low longitudinal 
momenta k Q <C k ~ 1/Zj_ , such that p c \/qf — is not negligible, and large positive 
scattering length < a ~ /j_ (meaning that a weak bound-state of V(r) approaches the 
threshold |Ï3|). For scattering by a central field, not only V(r) but also Í7(p) can be quite 
general. 

Viewing CIR as a low energy resonant scattering, one could say that bound-states close 
to threshold are neither probed at "high" energies k ~ l/l± (k — * q±, p c — > 0), nor do 
they exist for small scattering lengths (a d±). However, by calculating the bound- 
state with energy E' B , this interpretation for the physical mechanism behind CIR is not 
accurate: fo g ~ — 1 occurs before E' B approaches zero (threshold without confinement), 
whereas E' B — > eo (threshold under confinement) occurs only if l± is decreased much further 
below its CIR value. This is explicitly verified e.g. when U (p) is a square-well box of radius 
l± : using a cosine approximation to J for its roots, q ~ 37r/4/^ and Çi ~ 7n/4:l±, whence 
C = d±y/qj — q 1 = a/20/3 = 2.58 (see [ï(| for parabòlic U(p) and zero-range atom-atom 
interaction). 

In fact, the outer l = bound-state of V(r) in the absence of the confinement has the 

n 

energy Eb = — h 2 K B /2fi that is related to a via kb ~ 1/H, when a ^> Ry [14]. Under 
lateral confinement, its tail e~ KE>r is changed to be zero at the edge r = p = l±. By 
the uncertainty principle, this slight squeeze lifts Eb < by an amount eo, which can be 
sufficient for this state to pass the limit E = as l± decreases further. This new confmed 



bound-state E' B satisfies Eq.tJTJ) with k 2 replaced by 2p,E' B /h 2 , i.e., k = ±.i^Jq 2 — 2pE' B /h 2 . 
Since the diverging term e lk ° z should be absent from Eq. (|4ajl and e ífc °' 2 ' should decay, l//^ 
must vanish at Imjfco} > 0. From Eq. (jl3a)) . for a < 0, the virtual bound-state with energy 
Eb turns into a real one with energy E B , which starts at zero for a/d± = and goes to 
a positive fraction of eo as a/d± — > — oo. This bound-state exists only under confinement 
and its experimental measurement is reported in Q|. For a > 0, one obtains E' B —>■ Eb 
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for d± — > oc, as expected. For a — > +00 (or d± — > 0), tends to a positive fraction of 



e . It turns out that the CIR condition (at a/d± = 1/C = a/3/20 ~ 0.39) occurs before 
E' B reaches zero (at a/d± ~ 0.82). On the other hand, the CIR almost coincides with the 
condition E' B + (ei — e ) = e 

(a W ^ 0,5). In Re, Q, f* coincinence is eX ac t , 
since E 1 ^ + (ei — e ) can be associated with a bound-state of the excited channels n > 1 due 
to a special property of the harmònic oscillator. However, despite this coincidence, a general 
mechanism behind CIR needs further study, since E' B + (e x — e ) has no clear meaning yet 
beyond parabòlic guides and zero-range pseudopotentials. 

Excited Channels. At energies k 2 = q 2 E + k 2 lE > q 2 , the case is more complex. Keeping 
only the l = wave as before, the n-th scattering amplitude is 

ft - fn 9 f 7 b m N m /N n (i4) 



1 + a n + z 



d 2 

— ±(-kcotS - Poo) 



kr. 



where < m,n < n E , P 00 = (qf +nE - k 2 ) 1 ' 2 and in a n = Y^mjLn N mkn/N 2 k m , m = n is 
excluded. For the single incoming excited channel i.e., b n = 5 n ^ nE1 the amplitude f nE9 
does have the form Eq. (jl2|) at small k nE . Thus, CIR at threshold energies k — > q nE can occur 
when — tan5o/fc = d±/C as first indicated in Ref. [Uj for parabòlic confinement. In a more 
realistic situation of finite temperatures T, however, for a given energy each b n has the same 
weight (depending on E/T and with random phases). Since f ng = —b n cannot be met for 
all n simultaneously, one expects no sharp resonance, with the transmission and reflection 
probabilities being distributed among all channels according to Eq. ()10cj) . 

l-couplings. In the single mode regime, Eq. (|lüa|) is also an equation for t[ = T[/k 
without the singularity "jf ~ /cq . If then J^;/j;](2/' -|- l)Pii'ti> on the r.h.s converges, one can 
neglect it compared to ol\ for ko — > 0, and U pd ai/[i^kok — (21 + l)koPu ~ kok cot Si] is 
well behaved. Thus, angular momentum couplings should be negligible for ko — > and the 
series Eq. fll0b|) of individual momenta l is dominated by l = since 5i ~ k 2l+1 ~ 
are small, justifying Eq. ()13a|) . This does not apply straightforwardly to the excited channel 
case, whose approximation is based only on the smallness of Si . 

Discussion. Consider now the case U(p) = fiu)'j_p 2 /2 of harmònic confinement, /i being the 
reduced mass. In Eq. ()13bj) . the oscillator length a± = (h/ puj±Y^ 2 should replace d± = l±/N 
instead of Ix- This is due to tunneling, since \(p n (p)\ 2 ~ e ~p 2 l a \ i s small at p ~ l±_ (as 
in the square-well case) only if l± > a±. Then e\ — e = h 2 (qf — ql)/2p = 2hjj±_ and 
C = d±^qj — q 2 = 2. The difference to C — 1.4603... of Ref. |9j originates from the 



continuum limit in Eq.([5a|) and Eq.(|6bp. Indeed, from Eq.(9) of Ref. |9j], the continuum 
approximation for C is C = lim s ^ oo (J s ds'/y/s' — Yua'=i W 5 ') ~ lo^ 8 '/^ 8 ' = 2- In 
addition, this comparison reveals the nature of the "irregular" part 1/z of \l/(r) for the 
pseudopotential approximation (see Eq. (8) of Ref. [9j or the equivalent s-wave expansion in 
Eq.(9) of Ref. Jf|). This is the singular part of the free-space Green's function 7+Gi + 7-G2, 
with 7+ + 7_ = 1, and originates from the sum of the excited transverse levels. As a result, 
one expects certain details of the guide to be unimportant, except for the low lying levels 
which account for the terms 7 and A c and the bound-state E' B . 

We have provided a general treatment of quantum scattering in confined geometries. For 
scattering by obstacles inside the guide, the treatment should be applicable to a variety 
of central force fields V(r) and confining potentials U(p). For ultracold atòmic col·lisions, 
non-parabolic guides can be considered with restrictions due to the center of mass. The 1D 
scattering amplitude is given in terms of the free-space phase shifts ói and their couplings 
among each other. This covers the case of higher energies and a transversal multi-channel 
incident state. In the single mode regime, we have shown that the CIR is closely related to 
the behaviour of a confined bound state. 
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